
The Univarsity of Connecticut 
SCHOOL OF ENGINEERING 


Stow, Connecticut 06268 



AN EFFICIENT ALGORITHM FOR CALCULATION 


OF THE LUENBERGER CANONICAL FORM 


D. Jordan 
B. Srldhar 

Technical Report 72-9 



Department of Electrical Engineering 





AN EFFICIENT ALGORITHM FOR CALCULATION 


OF THE LUENBERGER CANONICAL FORM 

D. Jordan 
B. Sridhar 

Technical Report 72-9 


October, 1972 


This work has been sponsored in part by the 
National Aeronautics and Space Administration 
Research Grant NGL-07-002-002 



I ntroduc+ion 


The Luenberger canonical form is an extension of the control canonical 

[21 

form 1 - J for single- input or single-output controllable and observable systems 
to multivariable systems. The canonical form is not unique in the multivariable 
case. However, the controllability indicies are structural invariants of the 
system and correspond to the various blocks in the Luenberger canonical form.^ 3- ^ 
Consider the linear time-invariant controllable system 

x = Ax + Bu ( I ) 


where x is an nxl state vector and u_ ‘ Is an mxl input vector. In addition 


it is assumed that the columns of B are linearly independent. The controllability 

matrix r= [B,AB,A^B, . . . , A n *Bl has rank n and an nxn nonsingular matrix 
k I — 1 k -I k -I 

P = [b^| , Ab_|,...A b_|, b^ 2 » ^2'’'*^— 2 »•••> A_b m m 1 can be selected from 


the columns of T. Let e 

of P * where a. s .1 . k . . 

I J = l J 

transformat ion matrix 


, e 2 > • • .e m be the 0 | 
The vectors e| . 


. ..o^-th row respectively 

•e are used to construct the 
m 
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The transformation T reduces the system (I) to the form 

/■V A 

y = Ay + Bu 


( 2 ) 


where 
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There are several matrix computational difficulties, arising out of the 
need to find P * and T in arriving at the canonical form following Luen- 
bergers's construction. 

Here, a new algorithm Is presented which is more efficient and accurate 
than Luenberger’s construction. Also, the canonical form is computed directly 
and the transformation matrix T is computed only if necessary. 


Basic Resu I ts 


Let the transformation matrix 

T s Ct * ; t« . . . t • t . . . . t : ... t ] . 

_| _2 -c, -o, + 1, -<t 2 -o m 


(4) 


Then the similarity transformation satisfies the conditions 

AT = TA (5) 

B*TB (6) 

Numbering the columns of B as Cb_j bg ••• equation (6) imposes the following 

restrictions on the columns of T: 
b. = t 


b, = t + b^ , t 
-2 -o 2 a,, 2 -a, 


b = t + b t + . . . + b t 

— m — o o , , o t , m —0 . 

m I I m- 1 m- 1 

solving for , we can find constants such that 


—o 


* l C jl -J w 


i j=l 


j with c. . = I , i = I , . .m. 


(7) 


Further, from equation (5) the columns of T are related by the set of equations: 

m 

J = A \-J+i ", °a r a -J+l ±c, j = 1 , 2, . . . k - 1 ; £=l,2,...m (8.a) 

X/ 

m 


0 = At ... I 
— -a_-k.+ . L . 


W i=i 


1= 1 ,2, . . .m 


(8.b) 
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Examinlng (8) recursively. It can be seen that t^ and 0^ can be written entirely 
in terms of the t ’s by the equations 




. = A J 't 


m j-l 

l I a. 


A k t 


”°i 1 = 1 k=0 '°|'VJ +kf| -°1 


j = 1 ,2, . . .k -I ; £=l , . . ,m 


m V 1 


- I l 


A k t 


^£ k ^o V V k £ +k+l -I 


£=.1 , . . .m 


(9. a) 
<9.b) 


Substituting for t from equation (7) in equation (9.b) and rearranging the 


I 


terms, we get 


k £~* m m 


k £ a- in in i 

o=A C t c j£ ilp ~ I I I * a a a _ k + k +| • c j | * A 1-Lt £ =l * • • ,m 

j=l J J k=e 1=1 j = ? a J’°£ K £ K 1 ,J 

r si 

Kalman has suggested use of the decomposition form: 

. k.-l 

k m * u 

A b = - l l a A K b. , £= l,2,...m 


( 10 ) 


i e I k=0 


(II) 


where are constants and {A K b.} i = 1,2, ...m; k=l, ...kj is a basis for 

the n-space. Unfortunately, the accompanying change of basis suggested by 
l~5l 

Kalman 1- J does not lead either to the Luenberger canonical form or to the form 
suggested by him, as can be seen by applying his method to the example posed 
here! n. 

Without loss of generality we can assume that kj >_ ^ ^ ••• Then, 

k 

by the selection procedure of basis vectors A bj we can guarantee that 

a £ik = 0 for k > k £ a nd a £|k =0 for k =k £ , i > £. By defining a ££k =|, 

£ 

condition (II) can be rewritten as 

k -I 


0 = A 


l a o ik 
j=l * Jk £ J 


m £ 

1 

t=l k=0 


+ J £ ^ J“Lj * ^ m * (12) 
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Comparlson of equation (10) and (12) suggests the following procedure for 
selecting c's and a ? s: 


( I ) c = a„ .. 
J* *Jk 0 


a = 1,2,. . .m; j = I ,2, . . .£ 


" i+h C U " “Uk ( = 1 

m 

(,l) j=i C ij a c .a -k + k + I = " a £ik £ = 1,2, — m; I = 1,2, . . .m; 

J £- ^ 

■ ■ ■ k = 0,1,..,, k^- 1 

Given Cjj calculated In (I) and a^j k = 0 for k > k & , I ^ 

k = k £ , I > £ 

The a's can be obtained by backward substitution from (li). 

Also the b’s can be obtained In the same manner through the equation 

p - 1 

(III) b +2, c,. b + c. = 0, k = p-l, p-2, . . . 1 ; p=2,3..'.m. 

°k,p j=k+l kJ °j,p kp 

The algorithm suggested here computes the canonical form directly and, If 
necessary, the transformation matrix T can be obtained from equations (7) and 


(8. a) . 


Computational Efficiency 

The algorithm suggested in this paper results In a large reduction In the 
amount of computation necessary to obtain the canonical form. A comparison 
of Table I and Table 2 shows that, using Gaussian elimination techniques -where 
applicable, there is a saving of at least 2n^ multiplications (le about 43# 
reduction in the number of multiplications). 
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p 

n 2 ( n-m) 

-1 

4,3 , 

p 

y (n -n) 

T 

n 2 ( n-m) 

-| 

4,3 x 

T 

y (n -n) 

t“'at 

~ 3 
2n 

Ills, | 

nm 

Total 

4yn 3 + 2n 2 (n-m) - 4yn + nm 


o 2 3 •> 

! 

3 


Tatele I 

Total number of multiplications in 
Luenberger’s Method 


basis 

n 2 (n-m) 

a ! s 

y-~ y + m(n 2 -l ) 

a ij 

mn(m- 1 )/2 

b’s 

ro 

E 

1 

to 

E 

T 

3 3 2 2 1 

n + ^m -2n m + -^mn 

Total 

yTi 3 + ^ — 2mn(n-m) - y (4m+n) 


2 3 


< 2yn J 

Table 2 


Total number of multiplications In 
the new algorithm 


If only the canonical form is necessary there is a further reduction In the 
number of multiplications since the canonical form is computed directly without 
hav i ng to compute T . 

The algorithm should give better accuracy than Luenberger's method for 
two reasons. (I) the zeros are preserved in the canonical form and, as such, 
round-off errors in their computation are avoided. (2) the reduced amount of 
computation and Gaussian techniques lend to greater Inherent accuracy and the 
ability to refine the solution with additional computations. 


Examp I e 

Consider the system with 


'I 

0 

0 " 



‘l 0 ' 

0 

2 

0 

r 

B = 

0 1 

_0 

0 
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A= 
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Then, 



Thus kj = 2 and = I * Further, we have 




Th i s gives a | jq = 3, 

“ill " 4 ' 

01 120 

1 

a 2l 0 2' 

1 

a 2 1 1 ‘ ~V 

“220 


It Is assumed that oi| 2 |=0 and a 2\2 = ®’ 

Now, from ( I ) and ( i I \ Cj | = l , c \2~ a 2\\ = and c 22 = * * 

£ i k 

120 C 22 a a 2 ,o ) -k 1 + k + I “|20 * a 3l = 0 

110 c || a a ( ,a | -kj + k+ I +C I2 a o 2 ,a ( - k, + k + I * " °l 10 

0 2 1 2 * a 3 1 “3 I • g • 0 2 j “ -3 . 

121 c 22 a o 2 ,o ( -k, + k + I =_a 121 * a 32 = °* 
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C | I a O| t O|-k| + k + I + c |2 a cr2,0|-k| + k + 
* a 22 “ 2 a 32 = 4 1 ,e ’ a 22 = 4 


= -a. 


220 c OQ a . , . , 

22 a 2 >a 2~^2 + ^ + 


210 c. . a 


" a 220 a 33 = 2 


+ c 10 a 


II a | > a 2~^2 + ^ + I 12 °2’ a 2~^2 + k + 


9 23 * 2 a 33 = 0 1 ,e * a 23 = 2 


- -a 


210 . 


B can be obtained from (III). 
b_ o + c.„ = 0 I .e. b 


o, ,2 12 

This gives the canonical form 


0 | ,2 


= - c 


_ J_ 

12 2 ' 



' 0 

1 

0 ' 


' 0 

0 ■ 

A = 

-3 

4 

1 

2 

, B = 

1 

1 

2 


_ 0 

0 

2 


_ 0 

1 _ 


Application of (8. a) yields 


T = 


-3 I 
0 0 


'2 


Conclusions 

A new algorithm is suggested to obtain the Luenberger canonical form for 
multivariable systems. This method computes the canonical form directly with- 
out having to compute the transformation matrix. In addition, there Is a 
large reduction In the number of calculations. The reduced computations along 
with Gaussian techniques lend to greater inherent accuracy and the ability to 
refine the solution with additional computations. 
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